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1.0  Introduction 


This  document  is  the  final  technical  report  for  the  project  titled  "Development  of  Signal 
Processing  Algorithms  for  Space-Based  Radars”,  under  Contract  BAA-00-07-IFKPA,  which 
was  performed  by  Stiefvater  Consultants  for  the  Air  Force  Research  Laboratory,  Sensors 
Directorate  (SNRT),  located  in  Rome,  New  York.  The  objective  of  this  effort  was  to  develop 
signal-processing  algorithms  for  the  detection  of  ground  moving  targets  using  space-based 
radars  (SBR).  Particularly,  this  report  summarizes  our  work  on  the  along-track 
interferometric  synthetic  aperture  radar  (ATI-SAR)  approach  for  detecting  moving  targets. 


1.1  Background 

Historically,  Ground  Moving  Target  Indication  (GMTI)  has  been  accomplished  using 
manned  airborne  platforms,  such  as  the  Joint  Surveillance  Target  Attack  Radar  System  (Joint 
STARS).  Recently,  there  has  been  a  growing  interest  in  augmenting  these  airborne  assets 
with  a  space-based  capability.  Such  a  system  would  be  capable  of  providing  both  wide-area 
(theatre)  surveillance  and  tracking  of  airborne  and  ground  moving  targets.  The  space-based 
radar  (SBR)  capability  (Figure  1)  is  particularly  attractive  because  it  provides:  (1)  deep 
coverage  into  areas  typically  denied  airborne  assets;  (2)  greater  ease  and  flexibility  for 
deploying  the  sensor  platform  on  station  and  meeting  coverage  tasking;  (3)  greater  area 
coverage  rate  performance;  and  (4)  a  steep  lookdown  capability  for  foliage  penetration 
(FOPEN)  operation. 

Space-based  surveillance  requires  a  nominal  area  coverage  rate  of  several  hundred  km 7s 
with  revisit  rates  of  one  to  two  minutes,  while  tracking  requires  somewhat  shorter  revisit 
times.  Thus,  a  constellation  of  satellites  would  be  necessary  to  meet  these  requirements. 

Even  though  it  may  seem  that  the  altitude  of  a  satellite  can  be  freely  chosen,  the  two  Van 
Allen  radiation  belts  limit  practical  orbit  selection.  The  two  Van  Allen  radiation  belts  are 
centered  on  the  earth’s  geo-magnetic  axis,  at  altitudes  ranging  from  1500km  to  5000km  and 
from  13,000  to  20,000  km.  To  minimize  the  radiation  damage  to  electronic  components  in  a 
lightweight,  unshielded  satellite  design,  the  satellites  would  have  to  be  placed  in  orbits 
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outside  of  these  belts.  Therefore,  either  a  medium-earth  orbit  (MEO)  at  altitudes  of  .5000km 
to  13000km  or  a  low-earth  orbit  (LEO)  at  an  altitude  less  than  1500km  is  desirable. 

Today,  both  LEO  and  medium-earth  orbit  MEO  constellations  are  being  considered  for  SBR. 
Ultimately,  the  orbit  selection  will  depend  on  the  exact  system  requirements  and  the  total 
life-cycle  cost.  Although  smaller  numbers  of  satellites  would  be  required  at  MEO  altitudes 
than  at  LEO  altitudes  to  meet  the  large  surveillance  coverage  requirements,  the  low  life-cycle 
costs  (including  payload  plus  launcher)  make  LEO  satellites  more  attractive  for  regional 
coverage  than  their  MEO  counterparts.  This  effort  focuses  on  a  LEO  constellation.  The  high 
orbital  velocity  and  large  clutter  range  span  make  clutter  suppression  and  detection  of  slow- 
moving  targets  more  difficult. 


Figure  1.  Space-based  radar  concept:  satellite  coverage  provides  wide-area 
surveillance  and  tracking  of  airborne  and  ground  moving  targets. 


1.2  Space-Time  Adaptive  Processing 

The  ability  of  an  SBR  system  to  suppress  the  clutter  interference  is  complicated  by  the  large 
platform  velocities  associated  with  SBR  operation.  Such  operation  generally  requires  the 
detection  of  targets  within  the  clutter  Doppler  spectrum  (endo-clutter  case).  The 
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challenge  is  more  pronounced  at  LEO  deployments  where  orbit  speeds  are  much  larger.  For 
example,  the  mainbeam  clutter  velocity  spectrum  for  a  LEO  satellite  with  a  speed  of 
7,612  m/s  (at  an  altitude  of  500  km),  and  a  50-meter  antenna  in  azimuth  operating  at  L-band, 
is  about  ±36.5  m/sec.  In  other  words,  all  targets  with  speed  between  -36.5  to  ±36.5  m/sec  are 
immersed  in  main  beam  ground  clutter,  making  target  detection  using  conventional  pulse- 
Doppler  (PD)  processing  questionable. 

Three  decades  of  research  and  development  has  shown  that  space-time  adaptive  processing 
(STAP)  [1]  is  a  potentially  attractive  technique  for  clutter  suppression  in  airborne  radars.  It 
could  be  effective  in  SBR  applications  as  well  [2],  However,  space-based  STAP  is  a  much 
more  difficult  problem  than  its  airborne  counterpart  and  may  not  meet  the  requirement  on 
Minimal  Detectable  Velocity  (MDV),  because  of  the  high  platform  speed  and  large  clutter 
extension.  To  illustrate  this,  consider  an  L-band  radar  with  a  50  m  x  2  m  phased  array  carried 
by  a  low-earth  orbit  (LEO)  satellite  at  a  500  km  altitude.  Using  the  high-fidelity  Signal 
Modeling  and  Simulation  (SMS)  tool  [3],  we  simulated  the  clutter  covariance  matrix  with  the 
assumed  system  and  scenario  parameters  listed  in  Table  1. 


Table  1 

1.  System  and  scenario  parameters. 

Parameter 

Value 

Comments 

Antenna  Size 

50  m  x  2  m 

Number  of  Panels 

32 

No  overlap 

Panel  Size 

12  by  12  elements 

Center  Frequency 

1260  MHz 

Bandwidth 

10  MHz 

PRF 

2000  Hz 

Duty  Cycle 

8% 

Pulse  Number 

16 

Each  CPI 

Xmt  Weighting 

Uniform 

Rev  Weighting 

Uniform 

Uniform  within  subarray 

Platform  Fleight 

500  km 

Inertial  speed:  7612.61  m/s 

Orbit  Type 

Polar,  circular 

Elevation  Scan  Angle 

-28.3  0 

Mechanically 

Slant  Range  of  Interest 

1258.15  km 

Grazing  angle:  18.322° 
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Figure  2  shows  the  SINR  loss  (i.e.  loss  relative  to  the  clutter-free  case)  as  a  function  of 
Doppler  frequency  for  the  Matched  Filter  (MF)  [1]  and  the  joint-domain  localized  (JDL) 
algorithm  [4].  The  JDL  algorithm  is  a  post-Doppler  beam-domain  approach.  It  first 
transforms  the  radar  data  cube  into  the  beam-Doppler  domain,  and  then  performs  a  joint- 
domain  adaptation  over  a  N0  x  Jo  local  processing  region  (LPR),  where  N0  is  the  number  of 
Doppler  bins  and  Jo  is  the  number  of  spatial  beams.  The  low  degrees-of-freedom  (DOF) 
requirement  characterized  by  this  algorithms  leads  to  a  significant  reduction  in  required 
sample  support  and  computational  load.  Figure  2  shows  an  example  with  a  3x3  (3  Doppler 
bins  and  3  spatial  beams)  LPR.  The  corresponding  MDV,  measured  assuming  a  -5  dB  SINR 
loss  threshold,  is  13.81  m/s  for  MF  and  18.62  m/s  for  JDL.  For  a  space-based  GMTI  radar, 
the  requirement  on  MDV  could  be  under  lm/s.  Obviously,  a  STAP-based  approach  cannot 
meet  the  MDV  requirement  for  this  example  LEO  system.  Even  worse,  earth’s  rotation  will 
induce  a  crab  angle  to  the  platform,  which  makes  the  clutter  range-Doppler  spectrum  vary 
with  range  [5]  and  further  degrades  the  SINR  and  MDV  performance  of  STAP  [6]. 


Figure  2.  SINR  performance  for  MF  and  JDL. 
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1.3  ATI-SAR  Approach 

For  targets  moving  at  speeds  slower  than  the  endo-clutter  velocity,  alternative  GMTI 
techniques  may  be  needed  to  supplement  conventional  STAP  processing.  One  such 
technique  is  the  Along  Track  Interferometric  Synthetic  Aperture  Radar  (ATI-SAR*)  [7],  The 
ATI-SAR  technique  is  based  on  the  acquisition  of  two  complex  SAR  images  taken  under 
identical  geometries  separated  by  a  short  time  interval.  The  phase  difference  between  the  two 
interferometric  images  is  used  as  a  test  statistic  to  be  compared  with  a  decision  threshold. 

The  ATI-SAR  technique  has  been  proven  valuable  to  sense  the  earth-surface  motion  such  as 
ocean  surface  currents,  where  the  speed  accuracies  in  the  order  of  a  few  centimeters  per 
second  were  reported  from  airborne  platforms  [7].  Recently,  there  has  been  increasing 
interest  in  applying  ATI-SAR  techniques  for  the  diction  of  slow  moving  targets,  especially 
GMTI,  using  space-based  assets  [8].  A  major  objective  of  this  effort  is  to  assess  and  develop 
the  ATI-SAR  techniques  for  detecting  slow  moving  ground  targets. 

First,  the  detection  performance  of  conventional  ATI-SAR  is  examined.  It  is  shown  that  the 
high  false  alarm  rate  associated  with  this  technique  would  be  a  big  concern  for  any  military 
application.  To  reduce  the  false  alarm  rate,  a  dual-threshold  approach  that  combines  the 
conventional  interferometric  phase  detection  with  the  SAR  image  amplitude  detection  is 
proposed  in  this  work.  Pixels  with  strong  returns  could  contain  moving  targets,  stationary 
objects,  and  other  discretes.  The  amplitude  detection  suppresses  the  weak  returns  from  large 
smooth  surfaces  such  as  road  and  water  surfaces.  This  yields  two  important  results: 

(1)  An  interferometric  phase  map  (including  target  velocity  information)  obtained  by 
applying  the  interferometric  phase  detection  only  to  the  pixels  selected  by  the  amplitude 
detection, 


*  Also  known  as  AT  InSAR,  IFSAR  etc. 
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(2)  An  amplitude  map  (including  target  strength  infonnation)  obtained  by  applying  the 
amplitude  detection  only  to  the  pixels  selected  by  the  interferometric  phase  detection. 

We  have  illustrated  this  concept  [9]  using  the  Jet  Propulsion  Laboratory  (JPL)  L-band 
AirSAR  ATI  data,  collected  in  November  1998,  in  the  Monterey  Bay  area,  CA.  Further 
results  using  JPL’s  new  C-band  AirSAR  ATI  data  collected  in  April  2004,  in  the  area  of 
Lancaster,  CA,  are  presented  in  this  report. 

The  ATI-SAR  processing  also  requires  a  precise  calibration  of  the  platform’s  crab  angle. 
Because  the  data  were  collected  from  an  airborne  platform,  the  crab  angle  was  varying 
during  the  data  acquisition  flight.  Reference  [10]  developed  a  calibration  method  based  on 
the  inertial  navigation  unit  (INU)  measured  attitude  data  and  an  array  of  known  stationary 
comer  reflectors  (strong  scatterers)  as  the  reference.  In  this  report,  it  is  shown  that  the  crab 
angle  induces  a  range-dependent  modulation  to  the  interferometric  phase  and  a  blind 
calibration  method  that  does  not  require  any  knowledge  of  the  ground  reference  scatterer 
and/or  the  actual  crab  angle  is  introduced. 

The  new  C-band  AirSAR  ATI  data  we  received  are  collected  in  the  following  three 
modes[l  1],  which  will  be  explained  in  the  next  Section: 

•  Standard  mode:  Single  baseline,  single  transmitter. 

•  “Ping-pong  “  mode:  Zero  baseline. 

•  Double  baseline  mode:  Double  baseline. 

The  AirSAR  data  include  two  2.5  GB  image  data  files  and  many  supportive  files  for  each 
mode. 
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2.0  Basic  Measurement  Principles 


The  ATI-SAR  technique  is  based  on  the  acquisition  of  two  complex  SAR  images, 
taken  under  identical  geometries  separated  by  a  short  time  interval,  with  the  interferometric 
phase  being  used  as  a  test  statistic. 

Throughout  this  section,  we  assume  that  the  SAR  data  are  collected  in  an  ideal 
unsquinted  (zero-Doppler)  stripmap  mode  [12]  through  two  antennas  located  along  the  track 
with  separation  Bx  between  their  phase  centers.  We  further  assume  that  the  platform  moves 
along  a  straight  and  level  trajectory  at  a  height  Hp  and  speed  Vp.  An  approach  for  overcoming 
crab  angle  effects  will  be  presented  in  Section  5.0. 

Figure  3  shows  three  ATI  data  collection  modes  used  in  the  JPL’s  AirSAR  data.  The 
“standard”  mode  uses  a  single  transmitter  where  one  antenna  transmits  and  both  antennas 
receive.  The  “ping-pong”  mode  uses  dual  transmitters  where  each  antenna  alternately 
transmits  and  receives  its  own  echoes.  The  double  baseline  mode  also  uses  dual  transmitters, 
but  one  antenna  receives  echoes  transmitted  by  another  antenna  and  the  two  channels  of  data 
are  alternately  collected  on  adjacent  pulses. 

As  an  electromagnetic  wave  travels  a  round-trip  distance  2p  to  and  from  a  scatterer, 
the  phase  changes  by 

tp  =  — 4rcp/ A, ,  (1) 


where  X  is  the  wavelength  of  the  radiation.  The  Doppler  shift  of  a  scatterer  is 

f  _  1  dtp  _  2p 
D  2n  dt  X 


(2) 


where  p  is  the  range  rate. 

If  Image  B  is  taken  at  a  short  time  interval  At  lagged  to  Image  A  under  otherwise 
identical  conditions,  the  interferometric  phase  for  a  scatterer  can  be  expressed  as 

4n 

9  =  27rfDAt  =  — —  At  Vt  (3) 

X 


where  Vt  =  -p  is  the  radial  speed  of  the  scatterer. 
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With  a  threshold  rj  e ,  which  is  determined  by  the  required  false  alarm  rate,  the  test  statistic 
can  be  expressed  as 


[|®|  —  ’He*  Ho 

l|e|<r,9,  Hj 

And  the  MDV  can  be  expressed  as 

MDV  =  -  . 

4jiAt 


(4) 


(5) 


(c)  Double  baseline  mode:  B  receives  echoes  transmitted  by  A;  A 
receives  echoes  transmitted  by  B. 


Figure  3.  Three  ATI  data  collection  modes. 
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2.1  “Ping-Pong”  Mode 

In  the  “ping-pong”  mode.  Images  A  and  B  are  acquired  from  the  fore-antenna  and  the  aft- 
antenna,  respectively,  as  shown  in  Figure  4,  The  time  delay  is  defined  as  the  time  required 
for  the  aft-antenna’s  traveling  to  the  fore-antenna  position,  which  is  calculated  as 

At=Bx/Vp,  (6) 


where  Bx  is  the  baseline  distance  between  the  phase  centers  of  two  antennas,  and  Vp  is  the 
platform  velocity.  The  interferometric  phase  can  be  rewritten  as 


0  =  — ■ —Vt. 

X  v 


(V) 


The  maximum  unambiguous  ATI  velocity  can  be  found,  by  letting  0  =  2n ,  as 


V, 


unamb 


2B„ 


(8) 


The  ATI  system  can  measure  velocities  unambiguously  if  |Vt  |  <  Vunamb/2,andtheMDV 
becomes 


MDV  = 


>-Vpfin 

4kBx 


(9) 


VP 

Aft- Antenna  Fore- Antenna 


Fore- Antenna,  Image  A,  t  =  t0 

Aft- Antenna,  Image  B,  t  =  t0+At,  where  At  =  Bx  /  V 

4n  B 

The  interferometric  phase:  0  = - —  Vt ,  where  V,  is  the  target  radial  speed. 


Figure  4.  Principles  of  ATI-SAR  (“ping-pong”  mode). 
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2.2  Standard  Mode 

In  this  mode.  Image  B  is  the  same  as  that  obtained  in  the  “ping-pong”  mode,  but  Image  A  is 
formed  using  two  separate  antennas  (one  transmits,  another  receives). 

In  an  ideal  unsquinted  stripmap  mode  [12],  the  SAR  pixels  are  formed  from  zero-Doppler 
output.  From  the  geometry  in  Figure  5,  we  can  see  that  the  stripmap  SAR  image  using  two 
antennas  is  equivalent  to  that  obtained  with  a  single  antenna  located  in  the  middle  between 
two  antennas  (generally  at  the  cross  point  between  the  bisector  line  and  the  baseline). 
Therefore,  the  standard  mode  interferometric  SAR  is  equivalent  to  a  “ping-pong’  mode  with 
half  the  baseline  distance,  i.e.  Bx/2.  Thus,  the  interferometric  phase  can  be  rewritten  as 

0  =  — ■ —yt.  (10) 

XV* 


The  maximum  unambiguous  ATI  velocity  becomes 


V. 


unamb 


B. 


(ID 


Figure  5.  Equivalent  single  antenna  SAR  for  a  two-antenna  stripmap-SAR. 
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The  MDV  becomes 


MDV  = 


^Vpri9 

2ttBx 


(12) 


2.3  Double  Baseline  Mode 

The  double  baseline  mode  also  uses  dual  transmitters,  but  one  antenna  receives  echoes 
transmitted  by  another  antenna  and  the  two  channels  of  data  are  alternately  collected  on 
adjacent  pulses.  The  time  difference  between  two  images  is  one  pulse  repetition  interval 
(PRI),  i.e., 

At  =  T,  (13) 

where  T  =  1/PRF  is  PRI  as  PRF  is  the  pulse  repetition  frequency.  The  interferometric  phase 
can  be  rewritten  as 


47T 

0  =  — —  TV. . 
X 


(14) 


The  maximum  unambiguous  ATI  velocity  can  be  found,  by  letting  0  =  2 re ,  as 


V 


unamb 


X_ 

2T  ' 


(15) 


The  MDV  becomes 


MDV  =  ^ 
4kT 


(16) 
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2.4  Summary  and  Examples 

Table  2  lists  the  expressions  derived  above  for  three  ATI-SAR  modes.  Two  numerical 
examples  associated  with  JPL’s  AirSAR  system  [10]  are  presented. 


Table  2.  Summary  of  equations  for  three  ATI-SAR  modes. 


Mode 

Ping-Pong 

Standard 

Double  Baseline 

Interferometric  Phase  0 

x  v,  ’ 

x  vp  ’ 

4rt 

—  TVt 

X 

Unambiguous  ATI  Velocity 

X 

V unamb  (ui/ s) 

2BX 

Bx 

2T 

MDV  (m/s) 

^Vpri0 

^Vpp6 

Mle 

4ttBx 

2ttBx 

4kT 

2.4.1  L-Band  Example 

A  typical  set  of  parameters  associated  with  JPL’s  L-band  AirSAR  system  are: 

•  Vp  =  216m/s, 

•  X  =  0.2424  m, 

•  Bx  =  19.7736  m,  and 

•  PRF  =  420  Hz. 

Figure  6  shows  the  interferometric  phase  as  a  function  of  target  velocity  for  the  three  ATI 
modes  discussed  above.  Table  3  shows  the  corresponding  MDV  and  unambiguous  velocities. 
It  is  shown  that  the  ‘ping-pong’  mode  provides  the  best  MDV,  but  the  lowest  unambiguous 
velocity,  and  the  “double  baseline”  mode  provides  the  worst  MDV,  but  the  highest 
unambiguous  velocity.  It  may  be  necessary  to  combine  these  modes  and  probably  STAP- 
based  approach  to  cover  the  GMTI  velocity  range  of  interest. 


12 


^  ^  'U\ 


Figure  6.  Interferometric  phase  as  a  function  of  target  velocities  (L-Band). 


Table  3.  MDV  and  unambiguous  velocities  for  three  ATI-SAR  modes  (L-Band). 


Mode 

Ping-Pong 

Standard 

Double 

Baseline 

Unambiguous  ATI  Velocity  Vunamb  (m/s) 

2.9281 

5.8562 

15.9894 

MDV  (m/s)  for  r)9  =  1  radian 

0.2107 

0.4214 

8.1016 

MDV  (m/s)  for  r|e  =1.5  radian 

0.3161 

0.6321 

12.1524 

2.4.2  C-Band  Example 

A  typical  set  of  parameters  associated  with  JPL’s  C-band  AIRSAR  system  are: 

•  Vp  =  214.77  m/s, 

•  ^  =  5.67  cm, 

•  Bx  =  2.0794  m,  and 

•  PRF  =  564  Hz. 
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Figure  7  shows  the  interferometric  phase  as  a  function  of  target  velocity  for  the  three  ATI 
modes  discussed  above.  Table  4  shows  the  corresponding  MDV  and  unambiguous  velocities. 


Figure  7.  Interferometric  phase  as  a  function  of  target  velocities  (C-Band). 


able  4.  MDV  and  unambiguous  velocities  for  three  AT 

-SAR  modes  (C- 

3and). 

Mode 

Ping-Pong 

Standard 

Double 

Baseline 

Unambiguous  ATI  Velocity  Vunamb  (m/s) 

2.9281 

5.8562 

15.9894 

MDV  (m/s)  for  r)9  =  1  radian 

0.4660 

0.9320 

2.5448 

MDV  (m/s)  for  r|e  =1.5  radian 

0.6990 

1.3981 

3.8172 
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3.0  Detection  Performance  of  Conventional  ATI- 


SAR 


3.1  False  Alarm  Rate 

The  probability  of  false  alarm  (or  false  alarm  rate),  Pfa,  can  be  calculated  from  the 
probability  density  function  {pdf)  of  the  interferometric  phase  in  the  absence  of  target.  For 
simplicity,  consider  a  Gaussian  clutter  associated  with  a  homogeneous  background  and 
additive  white  Gaussian  thermal  noise.  In  the  absence  of  target,  the  interferogram  is  the 
product  of  two  complex  Gaussian  correlated  signals. 


3.1.1  Probability  Density  Function  of  Interferometric  Phase 

Assume  that  two  corresponding  pixels  in  Images  A  and  B,  xi  and  X2,  are  joint  circular 
Gaussian  variables  with  zero  mean.  The  joint  probability  density  function  (pdf)  is  given  by 

pdf  (w)  =  exp{-  w  HC  w}  (iy) 


where 


w  = 


X. 


The  correlation  matrix  can  be  written  as 
C  =  e{w 


}= 

fii 

yl  " 

j 

^2  j 

where 

11  =E[|Xl|2], 

12  =  E[x2  ],  and 

and  y  is  the  correlation  coefficient. 

The  interferogram  is  a  new  random  variable,  given  by 
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(18) 


(19) 


(20) 


(21) 


The  joint  pdf  of  the  magnitude  and  phase  of  v  can  be  found  to  be 


Pdf(|v|,4>)  = 


2  v 

TXl 2  1 

1- 

V 

y 

1 

exp 


2|y|v|cos(|) 


f 


K, 


v 


2|v| 

III- 


(22) 


where  K0(-)  is  the  modified  Bessel  function  of  the  second  kind  of  order  0. 

The  marginal  pdf  of  the  interferometric  phase  is  derived  by  Bamler  and  Hard  in  [14]: 


Pdf  (<|>)  = 


1- 


1 


2n  1  -  y  cos2  (j) 


1  + 


cos  c()  cos  1  (-  |y|  cos  4>) 


1-  y  “  cos2  (j) 


(23) 


which  is  fully  characterized  by  the  correlation  coefficient  |y|. 

In  the  presence  of  additive  white  noise,  the  above  pdf  formula  for  the  interferometric  phase 
holds  true  with  the  equivalent  correlation  coefficient 


y  = 


To 

1  +  1/CNR  . 


(24) 


where  yc  is  the  clutter  correlation  coefficient  and  CNR  is  the  clutter-to-noise  power  ratio. 


3.1.2  Numerical  Examples 

Figures  8,  10  and  12  show  the  pdf  of  the  interferometric  phase  in  the  absence  of  target,  for 
different  CNR  and  yc  =  0.98 , 0.99,  and  1,  respectively.  The  corresponding  false  alarm  rates 

vs.  the  phase  threshold  r|eare  shown  in  Figures  9,  11,  and  13,  respectively.  The  numerical 
results  are  listed  in  Tables  5  through  7.  For  example,  Pfa=  5.4x10'  foryc  =  1  and 
r)e  =1.5  radians,  which  is  a  very  high  false  alarm  rate  for  any  surveillance  radar. 
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Figure  8.  Phase  noise  pdf  for  different  CNR  values  with  yc  =  0.98. 


Figure  9.  False  alarm  rate  for  different  CNR  values  with  yc  =  0.98 
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Figure  10.  Phase  noise  pdf  for  different  CNR  values  with  yc  =  0.99. 


Figure  1 1 .  False  alarm  rate  for  different  CNR  values  with  yc  =  0.99. 
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Figure  12.  Phase  noise  pdf  for  different  CNR  values  with  yc  =  1 . 


Figure  13.  False  Alarm  probability  for  different  CNR  values  with  yc=  1 
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Table  5.  Pfa  vs.  CNR  and  phase  threshold  (yc  =  0.98). 


CNR\r|e  (rad.) 

0.5 

1 

1.5 

2 

2.5 

0  dB 

0.6737 

0.4312 

0.2724 

0.1654 

0.0855 

10  dB 

0.3025 

0.1173 

0.0594 

0.0327 

0.0162 

20  dB 

0.1095 

0.0339 

0.0162 

0.0088 

0.0043 

30  dB 

0.0804 

0.0241 

0.0115 

0.0062 

0.0030 

40  dB 

0.0773 

0.0231 

0.0110 

0.0059 

0.0029 

Table  6.  Pfa  vs.  CNR  and  phase  thresholc 

l(yc  =  0.99). 

CNR\rje  (rad.) 

0.5 

1 

1.5 

2 

2.5 

0  dB 

0.6712 

0.4281 

0.2698 

0.1636 

0.0846 

10  dB 

0.2852 

0.1082 

0.0545 

0.0299 

0.0148 

20  dB 

0.0763 

0.0227 

0.0108 

0.0058 

0.0029 

30  dB 

0.0441 

0.0127 

0.0060 

0.0032 

0.0016 

40  dB 

0.0407 

0.0117 

0.0055 

0.0030 

0.0015 

Table  7.  Pfa  vs.  CNR  and  ph 

ase  threshold  (yc  =  1). 

CNR\rje  (rad.) 

0.5 

1 

1.5 

2 

2.5 

0  dB 

0.668692 

0.424951 

0.267186 

0.161782 

0.083577 

10  dB 

0.266857 

0.099022 

0.049543 

0.027129 

0.013398 

20  dB 

0.039964 

0.011462 

0.005417 

0.002910 

0.001427 

30  dB 

0.004215 

0.001164 

0.000546 

0.000293 

0.000143 

40  dB 

0.000423 

0.000116 

0.000054 

0.000029 

0.000014 

3.2  Probability  of  Detection 

The  pdf  in  the  presence  of  target  is  required  for  computing  the  probability  of  detection  PD. 
Unfortunately;  this  pdf  is  not  analytically  available.  A  Monte  Carlo  simulation  was  used 
instead.  Figure  14  shows  the  pdf  plots  with  100,000  trials  for  three  different  values  of  signal- 
to-clutter  ratio  (SCR),  for  a  CNR  value  of  20dB  and  yc  =  1 ,  and  the  interferometric  phase 
induced  by  the  target  equal  to  2  radians.  The  corresponding  PD  plots  vs.  the  phase  threshold 
are  shown  in  Figure  15.  For  example,  if  r|0  =1.5  radians,  then  PD=  0.805  for  SCR  =  10  dB 
and  Pd  =  0.024  for  SCR  =  0  dB.  Therefore,  rather  high  SCR  is  required  for  reliable 
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detection.  Recall  that  (Table  7)  the  corresponding  false  alarm  rate  is  Pfa  =  5.4x10'  ,  which  is 
too  high  for  almost  any  military  radar  application. 


Figure  15.  Detection  probability  for  three  different  SCR  values  in  presence  of  a  moving 
target. 
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4.0  New  C-Band  AirSAR  ATI  Data 


4.1  Data  and  System  Parameters 

The  new  AirSAR  ATI  data  provided  by  JPL  recently  are  collected  at  C-band  in  April  2004, 
in  the  Lancaster  area,  CA,  as  shown  in  Figure  16.  The  data  we  received  included  three  file 
folders,  each  including  two  2.5GB  image  data  fdes  and  many  supportive  fdes,  as  listed  in 


Table  8. 

Table  8.  Information  of  received  C-band  AirSAR  ATI  data. 


File  Folder  Name 

OutCDA06 

OutCAB03 

OutCZA02 

Mode 

Ping-pong 

Standard 

Double  baseline 

Processing  Mode  Mnemonic 

CDA 

CAB 

CZA 

58  Mojave 
Moja\%o  ^.Airport 


Edwards. 


Edwards  AFB 


'Rosamond 


E  Avenue  J. 


E  Avenue  K 


Hillside 

Park 


Palrpdale 


Palmdale  Blvd 


Palmdi 


^Avenue 


Agua  Uuk» 
.Airpark! 


Pass  yjtlEmi 


Littlerock 
Recreation  Area 


Acton 


Lancaster,  California,  United  States 


0  km 


3 

/ 

North 

Edwards^  j 

Edwards  AF 
Auxilary 
North  Base 


Edwards  Air 
Force  Base 


E  Avenu 


Copyright  ©  1988-2003  Microsoft  Coro  and/or  *s  suppliers  All  richts  reserved  http  vMvav  mcrosoft  com/streets 

©  Copyright  2002  by  Geographic  Data  Technology.  Inc  All  nghts  reserved  ©2002  Navigation  Technologies  All  rights  reserved  Th»  data  indudes  information  taken  with  permission  from  Canadian  authorises©  1991-2002  Government  of  Canada  (Statistics 
Canada  and/or  GeomaScs  Canada),  al  nghts  reserved 


Figure  16.  Map  of  Lancaster,  CA. 
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Figure  17  shows  the  AirSAR  instrument  (panels  behind  wing),  mounted  aboard  a  modified 
NASA  DC-8  aircraft1.  Typical  system  and  scenario  parameters  associated  with  the  data  are: 

•  Platform  Velocity:  Vp  =  214.77  m 

•  Platform  Height:  Hp=  8.6934  km 

•  Average  Terrain  Height:  662  m 

•  Baseline  Distance:  2.0794  m 

•  PRF  =  546  Hz 

•  Range  at  the  First  Range  Bin:  Rmin=  8768.93  m 

•  Slant  Range  Spacing:  ARS  =  3.33 1  m 

•  Radar  Center  Frequency:  5.2875  GHz  (Wavelength:  X  =  5.67  cm) 

•  Chirp  Bandwidth:  40  MHz. 


Figure  17.  AirSAR  instrument  (panels  behind  wing)  mounted  aboard  a  modified 
NASA  DC-8  aircraft. 


1  More  details  can  be  found  at  http://airsar.jpl.nasa.gov/index_detail.html 
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4.1.1  “Ping-Pong”  Mode 

Each  SAR  data  file  has  a  total  of  167279  records.  Each  record  includes  2000  complex  (8 
byte  IEEE)  range  samples.  Records  are  divided  into  46  patches/blocks  separated  by  10  blank 
records.  The  sizes  (number  of  records)  and  the  first  record  indices  (“1st  Rec  #”)  of  46  patches 
are  listed  as  follows: 


Patch  # 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Size 

6546 

3910 

3908 

3910 

3910 

3909 

670 

810 

809 

809 

1st  Rec.  # 

6 

6562 

10482 

14400 

18320 

22240 

26159 

26839 

27659 

28478 

Patch  # 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

Size 

810 

3909 

3908 

3910 

3909 

3908 

3909 

3910 

3909 

3909 

1st  Rec.  # 

29297 

30117 

34036 

37954 

41874 

45793 

49711 

53630 

57550 

61469 

Patch  # 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

Size 

3909 

3910 

3909 

3909 

3910 

3909 

3909 

3908 

3910 

3909 

1st  Rec.  # 

65388 

69307 

73227 

77146 

81065 

84985 

88904 

92823 

96741 

100661 

Patch  # 

31 

32 

33 

34 

35 

36 

37 

38 

Size 

3909 

3910 

3910 

3910 

3908 

3909 

3910 

3908 

1st  Rec.  # 

104580 

108499 

112419 

116339 

120259 

124177 

128096 

132016 

Patch  # 

39 

40 

41 

42 

43 

44 

45 

46 

Size 

3909 

3908 

3908 

3908 

3909 

3908 

3908 

3908 

1st  Rec.  # 

135934 

139853 

143771 

147689 

151607 

155526 

159444 

163362 

4.1.2  Standard  and  Double  Baseline  Modes 

The  SAR  data  for  “standard”  and  “double  baseline”  modes  have  the  same  size  and  structure. 
Each  SAR  data  file  has  a  total  of  167272  records.  Each  record  includes  2000  complex  (8 
byte  IEEE)  range  samples.  Records  are  divided  into  46  patches/blocks  separated  by  10  blank 
records.  The  sizes  (number  of  records)  and  the  first  record  indices  of  46  patches  are  listed  as 
follows: 
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Patch  # 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Size 

6543 

3910 

3909 

3910 

3910 

3909 

670 

809 

810 

809 

1st  Rec.  # 

1 

6554 

10474 

14393 

18313 

22233 

26152 

26832 

27651 

28471 

Patch  # 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

Size 

810 

3909 

3908 

3909 

3910 

3908 

3909 

3910 

3909 

3909 

1st  Rec.  # 

29290 
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4.2  Stripmap  SAR  Images 

4.2.1  Cross-Range  Resolution 


The  cross-range  resolution  for  focused  aperture  stripmap  SAR  can  be  expressed  [13]  as 


PRF-RX 


2NVp 


(25) 


where  N  is  the  pulse  number  in  a  coherent  processing  interval  (CPI)  and  R  is  the  slant  range. 

Table  9  lists  some  results  for  PRF  =  546  Hz,  Vp  =  214.77  m,  and  X  =  5.67  cm.  It  is  usually 
desirable  to  maintain  a  square  resolution  in  stripmap  SAR  mode  [13].  Because  the  slant  range 
resolution  in  JPL’s  C-band  AirSAR  data  is  3.33m,  N  =  256  pulses  in  a  CPI  may  be  a  good 
choice  for  this  example  from  Table  9.  Therefore,  256  pulses  will  be  used  for  SAR  image 
processing  in  this  report  if  not  otherwise  specified. 
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Table  9.  Cross-range  resolution  (ARC)  in  meters. 


Slant  Range 

Rmin  =  8768.93  m 

Rmax  =  15430.99  m 

N  =  1024 

0.62 

1.09 

II 

CM 

K> 

1.23 

2.17 

N  =256 

2.47 

4.34 

4.2.2  SAR  Images 

Because  all  data  of  three  modes  are  recorded  in  the  same  scenario,  they  produce  similar  SAR 
images.  This  report  only  presents  the  results  from  one  data  file,  recorded  on  the  Antenna  1  in 
the  folder:  OutCDA06  (“ping-pong”  mode).  A  line-by-line  direct-convolution  method  is  used 
for  the  stripmap  SAR  imaging  with  N  =  256  pulses  in  a  CPI,  which  leads  to  a  square 
resolution.  Figure  18  shows  the  SAR  images  (in  dB  with  80dB  of  dynamic  range)  for  42 
patches.  Because  there  are  too  few  records  in  patches  7  through  11,  they  are  ignored  here. 

The  AirSAR  radar  starts  recording  Patch  1  from  the  area  of  Mojave  Airport,  CA,  and  then 
down  south  along  Route  14.  Comparing  the  SAR  images  from  Patch  1  to  Patch  13  with 
Figure  19,  newly  taken  bird-eye  picture  in  the  area  of  Mojave  Airport,  we  can  see  that  the 
SAR  images  match  the  picture  very  well. 
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Figure  18.  Stripmap  SAR  images  using  the  C-band  AirSAR  data.  Cont.) 
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Figure  18.  Stripmap  SAR  images  using  the  C-band  AirSAR  data.(Cont) 
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Figure  18.  Stripmap  SAR  images  using  the  C-band  AirSAR  data.(Cont) 
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Figure  18.  Stripmap  SAR  images  using  the  C-band  AirSAR  data.(Cont) 
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Figure  18.  Stripmap  SAR  images  using  the  C-band  AirSAR  data.(Cont) 
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Figure  18.  Stripmap  SAR  images  using  the  C-band  AirSAR  data. 
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Figure  19.  Newly  taken  bird-eye  picture  in  the  area  of  Mojave  Airport,  CA. 
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5.0  Blind  Calibration  for  Group  Phase  Shift  Induced 
by  Crab  Angle 


5.1  Group  Phase  Shift  Induced  by  Crab  Angle 

In  the  ideal  case,  two  ATI-antennas  are  aligned  with  the  moving  track.  However,  the 
crabangle  (yaw  and  pitch)  makes  one  antenna  offset  from  the  moving  track  of  the  other. 

Figure  20  illustrates  a  case  where  the  platform  has  a  yaw  angle  0C ,  which  leads  to  an  offset 
distance  in  cross-track  direction: 


Ay  =  Bxsin0c  , 


(26) 


where  Bx  is  the  baseline  distance.  The  cross-track  offset  of  two-antennas  will  induce  a  phase 
shift  between  Images  A  and  B,  called  group  phase  shift  in  this  report,  which  must  be 
calibrated  before  the  ATI-SAR  processing.  Because  the  platform’s  crab  angle  varies  during 
its  motion  [10],  the  group  phase  shift  will  vary  with  pixels  along  the  track.  Moreover,  we  can 
show  that  the  group  phase  shift  also  rapidly  varies  with  slant  range. 

From  the  scenario  in  Figure  20,  we  can  express  the  path  difference  from  a  ground  point  P 
(with  zero  Doppler)  to  two  antennas  as 

Ar  «  Ay  cos  cj)  =  Bx  sin  0C  cos  cj)  =  Bx  sin  0c  ^/l-(H/R)2 


where  cj)  is  the  grazing  angle  at  Point  P,  H  is  the  platform  height  from  the  local  earth  ground, 
and  R  is  the  slant  range  from  P  to  the  middle  of  two  antennas. 

In  the  “ping-pong”  mode,  the  group  phase  shift  due  to  the  path  difference  will  be 


Acp  =  ^ Ar  =  ^Bx  sin0c ^l-(WR)2 


for  “ping-pong”  mode. 


(28) 


Similarly,  the  group  phase  shifts  in  “standard”  and  “double  baseline”  modes  can  be  expressed 
as 


Acp  =  ~~Bx  sin  0C  -y/l  —  (H/R)2  , 

A 


for  “standard”  mode 


(29) 


and 
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(30) 


4  71  / - 

A(p  =  — VpTsin9<,-y/l-(H/R)2  ,  for  “double  baseline”  mode. 
k 

Using  the  parameters  given  in  subsection  4.1,  Figure  21  shows  the  group  phase  shift  as  a 
function  of  pixel’s  slant  range  with  yaw  angle  equal  to  5°  and  10°.  Figure  22  compares  the 
group  phase  shifts  of  three  modes  when  the  yaw  angle  is  equal  to  5°.  It  is  shown  that  the 
group  phase  shift  varies  rapidly  with  the  pixel’s  slant  range  for  all  three  modes.  The  phase 
variation  rate  with  range  is  the  fastest  in  “ping-pong”  mode,  and  the  slowest  in  “double 
baseline”  mode.  Therefore,  a  two-dimensional  (along  track  and  across  track)  interferometric 
phase  calibration  is  necessary  before  further  ATI-SAR  processing. 


Figure  20.  Across-track  offset  between  two  antennas  induced  by  the  crab 
(yaw)  angle. 
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Figure  2 1 .  Group  phase  shift  induced  by  yaw  angle  as  a  function  of  slant  range 
(“ping-pong”  mode). 


Figure  22.  Group  phase  shift  induced  as  a  function  of  slant  range 
(Yaw:  5°). 
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5.2  Two-Dimensional  Blind  Calibration  for  Interferometric 
Phase 

Reference  [10]  developed  a  calibration  method  based  on  the  INU-measured  attitude  data  and 
an  array  of  known  stationary  corner  reflectors  (strong  scatterers)  as  a  reference.  This  report 
proposes  a  blind  calibration  method  that  does  not  require  any  knowledge  of  the  attitude  and 
ground  reference  scatterers.  This  method  tries  to  estimate  the  group  phase  shift  as  a  function 
of  pixel’s  position,  and  then  calibrate  the  interferometric  phase  map. 

As  discussed  in  the  previous  section,  the  group  phase  shift  varies  with  the  slant  range  as  well 
as  the  cross-range.  It  is  difficult  to  estimate  the  two-dimensional  phase  shift  without  further 
information.  In  this  report,  we  assume  that  the  crab  angle  varies  slowly  with  time. 
Particularly,  we  assume  that  the  crab  angle  remains  the  same  in  adjacent  Ni  records,  which  is 
corresponding  to  the  platfonn  flight  distance:  NiVp/PRF.  For  example,  the  flight  distance  is 
about  80m  for  Ni  =  200,  Vp  =  215  m/s,  and  PRF  =  546  Hz.  Thus,  we  will  divide  the 
interferometric  phase  map  into  groups  first,  each  including  Ni  cross-range  pixels.  Then,  we 
estimate  and  calibrate  the  group  phase  shift  as  a  function  of  slant  range  for  each  group.  It  is 
found  that  one  estimation  and  calibration  is  not  enough,  because  the  estimation  of  group 
phase  shift  is  based  on  wrapped  interferometric  phases  (in  [-71,71] ),  and  the  compensated 
phases  are  wrapped  again  which  leads  to  a  new  group  phase  shift.  An  iterative  procedure  is 
used  in  this  report,  instead. 

Assume  that  the  size  of  the  interferometric  phase  map  to  be  calibrated  is  NSXNC,  where  Ns  is 
the  number  of  pixels  in  the  slant  range  dimension  and  Nc  is  the  number  of  pixels  in  the  cross¬ 
range  dimension.  The  proposed  calibration  method  is  described  in  the  following  steps: 

1.  Divide  the  interferometric  phase  map  into  [fix  (Nc/Ni)]  groups,  each  with  size  Ns*Ni 
except  for  the  last  one,  which  includes  the  remaining  (<  Ni)  cross-range  pixels. 

2.  For  each  group,  estimate  and  calibrate  the  group  phase  shift  as  a  function  of  slant 
range. 

There  are  many  possible  estimation  methods.  Our  proposed  method  is 
described  in  the  following: 
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1.  For  each  slant  range  cell,  sort  the  interferometric  phases  of  Ni  cross¬ 
range  pixels 

2.  Take  the  mean  of  10  middle  phases  as  the  estimation  of  the  group 
phase  shift  for  the  corresponding  slant  range  cell. 

3.  Compute  the  mean  and  standard  deviation  of  the  estimated  group 
phase  shift  (over  slant  range  cells).  The  mean  and  standard  deviation  in 
the  ideal  case  are  zero.  If  they  are  small  enough,  no  further  calibration 
is  necessary,  and  go  to  Step  3. 

4.  Compensate  the  phase  shift  by  subtracting  the  estimated  group  phase 
shift,  and  wrap  the  results  within  [-71,  ti]  . 

5.  Back  to  Step  1). 

6.  Back  to  Step  2  for  next  group. 

5.3  In  Illustration  Example 

To  illustrate  the  problem  of  group  phase  shift  and  the  proposed  calibration  method,  we  have 
selected  a  partial  area  (500  slant  range  cells)  from  Patch  3  in  “OutCDA06”,  as  shown  in 
Figure  23  for  the  SAR  image  of  selected  area. 

Figure  24  shows  the  interferometric  phase  map  for  the  SAR  image  shown  in  Figure  23 
before  phase  calibration.  Obviously,  the  phases  of  most  pixels  are  far  from  zero,  as  shown  in 
Figure  25  for  the  histogram  of  interferometric  phase  before  phase  calibration. 

Figure  26  shows  the  interferometric  phase  image  after  the  phase  calibration  using  the 
proposed  blind  calibration  method.  Clearly,  the  phases  of  most  pixels  are  close  to  zero  after 
the  calibration,  as  shown  in  Figure  27  for  the  corresponding  histogram  of  interferometric 
phase  after  the  calibration. 
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Across  Track  Range  (km) 


Figure  23.  Image  A  for  selected  area  (from  Patch  3). 


Figure  24.  Interferometric  phase  map  before  the  phase  calibration  for  crab  angle. 
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Figure  25.  Flistogram  of  interferometric  phase  before  phase  calibration. 


angle(AB*)  After  Calibration 


Across  Track  Range  (km) 


Figure  26.  Interferometric  phase  map  after  the  phase  calibration. 
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Figure  27.  Histogram  of  interferometric  phase  after  the  phase  calibration. 
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6.0  Moving  Target  Detection  Using  ATI-SAR 


In  this  section,  we  still  use  the  selected  partial  area,  as  shown  in  Figure  23,  from  Patch  3  in 
OutCDA06  (“ping-pong”  mode),  except  where  otherwise  specified. 

6.1  Phase-Only  Detection 

Phase-only  detection  is  the  conventional  ATI-SAR  method.  Denote  r|e  as  the  phase 
threshold.  The  phase  at  any  pixel  below  the  threshold  will  be  forced  to  zero.  Figure  28  shows 
the  results  with  r|0  =  1  radian.  Obviously,  the  false  alarm  rate  is  too  high  for  radar 
applications  involving  Ground  Moving  Target  Indication  (GMTI). 


Figure  28.  Interferometric  phase  map  with  r)0  =  1  radian. 

6.2  Amplitude-Only  Detection 

The  amplitude-only  detection  suppresses  the  weak  pixels  from  large  smooth  surfaces  such  as 
road  and  water  surfaces.  This  detection  is  similar  to  the  conventional  constant  false  alarm 
rate  (CFAR)  processing.  There  are  many  algorithms  available  in  the  literature  to  determine 
the  threshold  [15].  The  performance  depends  on  the  environment.  As  an  example,  the 
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threshold  is  counted  relative  to  the  mean  amplitude  for  background  in  this  report,  which  is 
estimated  by  taking  the  root  mean  square  (rms)  mean  value  of  the  median  amplitudes  (of  the 
corresponding  slant-range  pixels)  over  cross-range  pixels. 

Figure  29  shows  the  SAR  amplitude  image  trimmed  with  an  amplitude  threshold  (r)a  =  6dB  ) 

only.  This  image  shows  stronger  pixels  that  could  possibly  contain  moving  targets,  stationary 
objects,  and  other  discretes.  The  amplitude-only  detection  suppresses  the  weak  pixels  such  as 
those  corresponding  to  road  and  water  surfaces,  but  it  can  not  separate  the  moving  targets 
from  the  stationary  background. 
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Figure  29.  SAR  amplitude  image  with  r\d  =  6  dB  . 


6.3  Dual-Threshold  Detection 

In  this  effort,  we  propose  to  combine  the  above-mentioned  amplitude-only  detection  with  the 
phase-only  detection.  Using  the  two  thresholds  (phase  and  amplitude),  we  get  two  outputs: 

An  interferometric  phase  map  (target  velocity)  is  obtained  by  applying  the  interferometric 
phase  detection  only  to  the  pixels  selected  by  the  amplitude  detection.  Specifically,  the  phase 
at  a  pixel  will  be  forced  to  zero  if  its  image  amplitude  is  below  a  pre-determined  threshold. 
And, 
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An  amplitude  map  (target  strength)  is  obtained  by  applying  the  amplitude  detection  only  to 
the  pixels  selected  by  the  interferometric  phase  detection. 

Figure  30  shows  the  results  of  the  interferometric  phase  map  obtained  by  applying  the 
amplitude  detection  results  of  Figure  29  onto  those  of  Figure  28,  i.e.,  forcing  the  phase  at  a 
pixel  to  zero  if  its  amplitude  is  less  than  the  mean  value  by  an  amplitude  threshold 
(r)a  =  6  dB  here).  Clearly,  the  false  alarm  rate  is  dramatically  reduced  using  this  approach. 

Similarly,  we  can  obtain  an  amplitude  map  of  potential  moving  targets  by  applying  the  phase 
detection  results  of  Figure  28  onto  those  of  Figure  29,  i.e.,  eliminating  the  pixels  whose 
phases  are  below  the  phase  threshold  (r|e  =1  radian  here),  as  shown  in  Figure  31,  which 

corresponds  to  the  interferometric  phase  map  of  Figure  30.  In  other  words,  Figures  30  and  31 
show  velocity  and  strength  information  of  potential  moving  targets.  The  locations  shown  in 
these  Figures  are  shifted  from  their  real  ones  due  to  the  SAR  processing.  It  is  possible  to 
restore  the  real  locations  of  those  slow  targets  without  any  interferometric  phase  ambiguity. 
Figure  32  further  combines  the  detection  results  with  the  original  SAR  image.  The  red  points 
along  the  road/highways  are  potential  targets  (moving  vehicles).  The  red  points  may  have 
shifted  from  the  road  tracks  due  to  the  fact  that  a  moving  target  appears  in  a  SAR  image  at  an 
apparent  position  shifted  in  the  along-track  direction  from  its  true  position. 

The  road  is  almost  in  parallel  with  the  flight  track  of  the  radar  platform.  Therefore,  the 
moving  vehicles  on  that  road  should  have  very  low  radial  speed.  The  fact  that  many  targets 
on  that  road  are  detected  shows  the  extraordinary  detection  capability  of  the  ATI-SAR 
approach  to  slow  moving  targets. 

A  MATLAB  file,  DTATRSAR.m,  is  given  in  Appendix  3.1  for  demonstrating  the  results 
shown  in  this  section. 
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Figure  30.  Interferometric  phase  map  with  r)e  =  1  radian  and  pa  =  6  dB  . 
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Across  Track  Range  (km) 


Figure  32.  Positions  of  potential  targets  (red  points)  on  the  SAR  image 
(“ping-pong”). 


6.4  Applying  the  Dual-Threshold  Approach  to  Other  ATI- 
Mode  Data 

The  above  demonstration  example  is  based  on  the  “ping-pong”  mode  data.  In  this  subsection, 
we  will  repeat  the  above  study  using  the  “standard”  and  “double-baseline”  mode  data,  taken 
from  file  folders  “OutCAB03”  and  “OutCZA02”,  respectively,  in  the  same  area  as  shown  in 
Figure  23. 

6.4.1  “Standard”  Mode  Data 

This  mode  is  similar  to  the  “ping-pong”  mode,  but  with  higher  MDV  as  shown  in  Table  4. 
Figure  33  shows  the  interferometric  phase  map  after  phase-only  detection  with 

r|e  =  1  radian.  Figure  34  shows  the  SAR  amplitude  image  trimmed  with  an  amplitude 
threshold  ( pa  =  6  dB )  only.  Figures  35  and  36  show  the  interferometric  phase  map  and  the 

amplitude  map  of  detected  targets  using  joint  amplitude  and  phase  detection,  respectively. 
Figure  37  shows  the  detection  results  on  the  original  SAR  image.  Compared  with  Figure  32, 
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we  can  see  similar  results  but  several  slower  targets  are  missed  because  of  the  higher  MDV 
in  the  “standard”  mode. 
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Figure  33.  Interferometric  phase  map  with  r|9  =  1  radian  (“standard”). 
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Figure  34.  SAR  amplitude  image  with  r)a  =  6  dB  (“standard”). 
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Figure  35.  Interferometric  phase  map  with  p0  =  1  radian  and  pa  =  6  dB  (“standard”). 


Figure  36.  Amplitude  map  with  r|e  =  1  radian  and  r)a  =  6  dB  (“standard”). 
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Figure  37.  Positions  of  potential  targets  (red  points)  on  the  SAR  image  (“standard”). 


6.4.2  “Double  Baseline”  Mode 

This  mode  has  the  worst  MDV  as  shown  in  Table  4.  Figure  38  shows  the  interferometric 
phase  map  after  phase-only  detection  with  r)e  =  1  radian.  Figure  39  shows  the  SAR 

amplitude  image  trimmed  with  an  amplitude  threshold  (r|a  =  6  dB )  only.  It  is  shown  that 

both  phase-only  and  amplitude-only  detections  produce  a  very  high  false  alarm  rate,  just  as 
those  seen  in  the  “ping-pong”  mode  and  “standard”  mode.  Figures  40  and  41  show  the 
interferometric  phase  map  and  the  amplitude  map  of  detected  targets,  respectively,  using 
joint  amplitude  and  phase  detection.  Figure  42  shows  the  detection  results  on  the  original 
SAR  image.  Compared  with  Figure  32,  we  can  see  that  most  targets,  detected  in  the  “ping- 
pong”  and  “standard”  modes,  have  disappeared  because  the  “double  baseline”  mode  can  only 
detect  high  radial  speed  targets. 
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Figure  38.  Interferometric  phase  map  with  r)e  =  1  radian  (“double  baseline”). 
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Figure  39.  Amplitude  image  with  r)a  =  6  dB  (“double  baseline”). 
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Figure  40.  Interferometric  phase  map  with  r|9  =  1  radian  and  r)a  =  6  dB 
(“double  baseline”). 


Figure  41.  SAR  Amplitude  map  with  r|e  =  1  radian  and  r|a  =  6  dB 
(“double  baseline”). 
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Across  Track  Range  (km) 


Figure  42.  Positions  of  potential  targets  (red  points)  on  the  SAR  image 
(“double  baseline”). 
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7.0  Conclusions 


Conventional  ATI-SAR  approaches  can  detect  targets  with  very  low  radial  speeds,  but  their 
false  alann  rate  is  too  high  to  be  used  in  GMTI  radars.  The  proposed  dual-threshold 
approach,  which  combines  the  conventional  interferometric  phase  detection  and  the  SAR 
image  amplitude  detection,  can  effectively  reduce  the  false  alarm  rate.  The  concept  of  the 
dual-threshold  approach  is  illustrated  using  JPL's  AirSAR  ATI  data.  A  simple  two- 
dimensional  blind-calibration  procedure  is  proposed  to  correct  the  group  phase  shift  induced 
by  the  platform’s  crab  angle.  However,  the  work  presented  in  this  report  is  only  our  very 
early  effort  in  this  area.  Future  work  would  include:  (1)  completing  the  dual-threshold 
detection  theory  with  performance  analysis  that  leads  to  an  easy  determination  of  thresholds; 
(2)  further  work  with  the  AirSAR  ATI  data,  including  the  completion  of  blind-calibration 
method,  comparative  studies  in  different  modes  with  more  sample  images,  especially  those 
images  with  known  targets;  (3)  assessing  the  SAR-MTI  algorithm  suggested  in  [16]  using  the 
AirSAR  data;  (4)  modifying  the  dual-threshold  approach  by  replacing  the  SAR  amplitude 
detection  with  the  SAR-MTI  detection. 
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